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Abstract 

The equivalence between Tsallis Thermodynamics and Hill's Nanothermodynamics 
is established. The correct thermodynamic forces in Tsallis thermodynamics are es- 
tablished. Through this connection we also find a general expression for the entropic 
index q which we illustrate with two physical examples, allowing in both cases to 
relate q to the underlying dynamics of the Hamiltonian systems. 
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1 Introduction 



During the last decade there has been growing interest in Tsallis (nonexten- 

sive) thermostatistics [l[ 0] . This formalism has been successfully applied to a 
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wide variety of statistical systems at scales ranging from particle and nuclear 



)liysics p| to astrophysics 



^ and situations involving long range interactions 



5| as well as low dimensional maps and multifractals 6]. Another important 
success is the straightforward explanation of the occurrence of Levy distribu- 



tions in nature 



7|. The building block of Tsallis thermostatistics, known as 



Tsallis entropy, was introduced in 1988 1| 

- 1 

5* = (g G 3?) 



l-q 



Here pj is the probability of the microstate j and S* is a generalized entropic 
measure which reduces to the Gibbs-Shannon one S = —k Y,^ Pj Inpj when the 
entropic parameter q tends to unity. The main feature of this entropic form is 
its nonadditivity. Given two subsystems _A and 5 of a composite nonextensive 
system A + B it can be easily checked J^] that 



(2) 



S*{A + B)= S*{A) + S*{B) + (1 - q)S*{A)S*{B) 



Despite the wide variety of applications that the formalism has proven 
there are some important points concerning Hamiltonian systems that have 

not yet been elucidated. Among these, the problem of how can q be calculated 

P 

for these systems has attracted great attention 8] and is, arguably, the main 
opened issue of Tsallis thermodynamics (TT) 9]. Another point that has not 
arified is what is the correct form of the thermodynamic forces in TT 



Deen c 



10 



11 



12| . It has also been stated that Tsallis entr opy cannot have a well 



defined thermodynamics because of its nonadditivity 



IG 



13|. 



Quite interestingly also, there exist already another previous important gener- 
alization of traditional thermostatistics which is well grounded physically and 

I I 

which was originally conceived for dealing with nonextensive systems [14] . This 



2 



formalism is Hill's Nanothermodynamics 



15| (NT) and was originated in 1962 



generalizing the concept of Gibbs' chemical potential to complex nonextensive 



systems that can vary its entropy by fragmentating in smaf 
las been applied succesfully since then to ferromagnetism 



3 



17| and liquid-vapor interfaces 



m 



er subsystems. NT 



16|, glassy systems 



19|. NT is well grounded physically from 



its very beginning in contrast with TT which relies in the ad hoc (although 
beautiful) postulate for the entropy inspired in the multifractal formalism 1|. 
TT has, however, many practical advantages and allows to deal quite straigth- 
forwardly with complex systems (specially those exhibiting power laws statis- 
tically since these can be closely described by means of Tsallis q-exponential 
distributions). NT could be misunderstood as applicable only to small sizes or 
small numbers of particles. But stars, clusters and multifractals can be consid- 
ered also as NT systems (the range of the interactions and/or correlations are 
of the order of the size of the system or longer) and, in general, all systems in 
which their equilibrium properties depart from the standard description, which 
relies in the concept of extensivity . NT is a rigorous statistical foundation 
of thermodynamic finite-size effects as well as others coming from the com- 
plexity of the interactions. It is therefore interesting questioning whether if 
there is any connection between TT and NT. 



In this article we clarify all these points. Specifically: i) We establish the equiv- 
alence between TT and NT, ii) this connection allows us then to introduce 
an adequate form for the thermodynamic forces in TT, overcoming previous 
difficulties 1^ iii) the nonadditivity property is then understood for complex 
systems thermodynamically and, importantly, iv) we provide an expression 
for q that can be used to calculate this quantity in Hamiltonian many particle 
systems relating it to its internal dynamics. 



3 



2 Formal equations of Tsallis Thermodynamics (TT) 



We present first previous rigorous results of TT and NT to better clarify 
our approacli. Maximization of the Tsallis entropy imposing the constraints 
{Xa)q = Y.f Pj^aii) / Yl^ p] fo^ the biased average {Xa)q of an extensive quan- 
tity Xa leads to recover a Gibbs-like differential equation for the entropy 

dS* = Y.y^d{X^)^, (3) 

a 

The rest of the formal equations of TT have been recently established 
The Euler equation of TT is 

= ^ [1 + (1 - q)S*] In [1 + (1 - q)S*] . (4) 

and the Gibbs-Duhem equation has the form 

Y,{X^)qdyo. = In [1 + (1 - q)S*] dS* (5) 

a 

A reasonable condition (local additivity) that variables {Xa)g must satisfy and 
which is used in the derivation of Eqs. (4) and (5) is d{Xa{XA))q/dX\^^^ = 
{Xa{A))q for a system A times bigger than another A. Note that this by 



111 ], since 



no means imply that {Xa{XA))q = X{Xa{A))q (global additivity) 
A is replaced by unity. Global additivity generally holds only when an en- 
semble of A independent noninteracting systems is considered (see below) 
but not necessarily upon subdivision of a system into smaller fractions (it 
is necessary, for example, to provide some energy to separate two droplets 
A and B from a bigger droplet A + B because of the energetic cost of cre- 
ating additional interfaces). Tsallis entropy satisfies JJ,J dS*{\A)/d\\y^^^ = 
[1 + (1 - q)S*{A)] In [1 + (1 - q)S*{A)] /(I - q). It is to be noted that, as oc- 
curs with their respective counterparts in standard thermodynamics, Eqs. (3) 



4 



to (5) are not independent, i.e. knowing two of them the other one can be 
trivially deduced. The above general equations and the Legendre transform 
mechanism allows to determine all thermodynamic quantities of a nonexten- 
sive system. 



3 Hill Nanothermodynamics (NT) 



Let us now introduce Hill's formalism of NT. It is based in maximization of 
the Gibbs-Shannon entropy but contains a new entropic thermodynamic po- 
tential, the subdivision entropic potential (SEP) J ^ which couples nonlinearly 
the natural thermodynamic variables of the system. Introduction of this new 
potential is made at the ensemble level, in which a set of M identical, noninter- 



acting small systems is considered 15|. In this case, the set of formal equations 
are J dM ^Y^aV^M {^c^u^t = ^St (Gibbs); JM + Y.aya,H {Xa)H,t = (Eu- 
ler) and NdJ + J2a i^a) h t^ya,H = (Gibbs-Duhem). Here, subindex H 
means "Hill's variables" which are the physical (averaged) locally extensive 
{{Xa)m) and intensive {ya,H) ones, t means "total" and denotes the prop- 
erties of the whole ensemble. Thus St is the total entropy for the ensemble 
of N systems. We have St = J\fS and {Xa) m = N {Xa) ^ in terms of the 
respective variables for one system (note that here the members of the ensem- 
ble do not interact directly although, of course, can exchange heat, volume 
and particles depending on the ensemble considered, this being itself a formal 
construction in which total properties are M times those of one system). The 



formal equations of NT become then for one system 



5 



J2ya,Hd {Xa)H = dS (6) 

a 

Y.y-,H{x^)H = s-J (7) 

a 

J2iXc.)Hdya,H = -dJ (8) 
a 

4 Connection between TT and NT 

Eqs. (6) to (8) are to be compared with Eqs. (3) to (5) respectively. It can be 
seen that if we make the following connection S = S*, ya,H = Va, {Xa) ^ = 
(X^)g and 

j.S--|l + (l-.)5l!^^^^^^ (9) 

the structure of the formal equations of TT and NT is the same. Here the 
Tsallis entropy S* is not only the physical one: its property of nonadditivity, 
see Eq. (2), is also the basis for the SEP J' which is found to be necessary 
to explain the thermal behavior of small systems (at least). TT describes, 
thus, the most general thermal equilibrium, the nanot her mo dynamic equilib- 



rium j20[, in which the new potential plays a decisive role. vanishes 
for a macroscopic system, for which one has also g = 1 in Eq.(9), and is a 
measure of the (thermodynamic) smallness of the system. In making this cor- 
respondence it could seem counterintuitive equating an additive entropy S to 
a nonadditive one S*. However, 5* is Gibbs-Shannon like because Hill's entropy 
is written statistically in terms of thermal averages considering a nanosystem 
in equilibrium with a macroscopic heat bath at a given average temperature. 
The nonlinear coupling of the intensive variables is accounted for thermody- 
namically in Eq.(7) by means of the generalized potential J'. This is totally 
equivalent to consider the nanosystem at the same average inverse tempera- 
ture arising this time not from coupling to a macroscopic heat bath but to a 



6 



neighbouring environment with fluctuating temperature following a chi-square 
distribution (without need of introducing J in this case). The latter leads to 



Tsallis statistics 



2l|, the entropy being nonextensive and the probability 
distributions depending on the parameter q (that comes from the order of the 
chi-square (Gamma) distribution considered). Some authors have shown that 
;he Tsallis distribution is the Laplace transform of the chi-square distribution 



8 



2l|. This is the key point in passing from Boltzmann-Gibbs distributions 



to Tsallis q-exponential distributions. 

Nonadditivity in the entropy (a statistical approach) and introduction of the 
potential J (a thermodynamic approach) are only different pathways that 
lead to the same thermodynamics. As shown below, thermodynamic excess 
functions (Hill's NT) are linked to entropic nonadditivity (Tsallis TT). It 
is also to be noted that the ensemble approach followed by Hill to derive the 
thermodynamics of a small system is equivalent to assume that local additivity 
of parameters and Eq. (6) hold as well as the following expression for the 
additivity breaking of the entropy dS{\A) j d\\-^^^ = S{A) — J . 

Despite the different statistical behavior both Hill and Tsallis entropies can be 
made to coincide numerically and the differences of both formalisms concern 
only the way of averaging. It is worthy insisting on the fact that thermodynam- 
ically speaking, however, these differences are irrelevant, since the structure of 
the formal equations is the same. Both formalisms, TT and NT, preserve the 




Gibbsian structure of thermodynamics [l|,ll5[. Fluctuations are also gaussian- 
like in both formalisms jlll-lisl. Furthermore, if the range < g < 1 is consid- 



ered to be the physically meaning 
g 7^ 1 than in the extensive case 



ul one fluctuations in TT are larger for 
ll| which is again consistent with the fact 



that fluctuations in NT for a system with 7^ are also larger than in the 



7 



J = Q case [15']. Although Tsalhs probabihty distributions are more general 
than Boltzmann-Gibbs ones, it is important to note that, in NT, J intro- 
duces additional degrees of freedom in the evaluation of the relevant partition 
function 16| (which in NT is always related to excess quantities besides the 
traditional thermodynamic ones). Furthermore, and quite importantly, both 
TT and NT share the property of ensemble non-equivalence for q ^ 1 and 
JT" 7^ 0, respectively. 



5 Correct form for the Lagrange parameters 



It has been questioned the physical meaning of the parameters t/q, and some 



authors 



11 



12| have concluded previously that the physical quantities of inter- 
est are ija = |/q,/[1 + (1 — This conclusion, however, is wrong. Parameters 
ija are naturally related to a thermodynamic formalism in which use of Renyi 
entropy (and not Tsallis entropy) is made. As a function of Tsallis entropy, the 



Renyi entropy S* has the form 5* = In [1 + (1 — 



while Tsallis entropy is nonextensive and stable 



but thermodynamically unstable 



standard thermodynamics 



a)S*] /(I - q) 



l|. However, 



22j | , Renyi entropy is extensive 



and cannot be employed to generalize 



22| . Parameters ya also arise when it is assumed ad 



hoc global additivity of energy and the maximum Tsallis entropy for two sys- 



tems brought into contact is considered 



12| . This, however, is also incorrect. 



since energy cannot be globally additive where entropy does not satisfy this 
property (otherwise, temperature could not be an intensive variable 10|). It 
is to be noted therefore that neither the correct thermodynamic forces nor the 
physical meaning of the formal relations were previously established. Through 
the above correspondence and from Hill's NT, it is now known that the Ua 



8 



must be equal for two different systems put in contact at equilibrium. The 
Ua are then the physically meaningful thermodynamic forces. This allows to 
establish equilibrium TT, which is now free from recent criticisms ^ arising 
from the use of incorrect thermodynamic forces. 



6 Excess quantities. Meaning of the nonadditivity of Tsallis en- 
tropy 

Excess thermodynamic quantities {Xa)'f \ as considered, for example, in sur- 



face physics 



15 



2|, can be defined as (X«(A,5))J^) = + B))^ - 

{Xa{A))q — {Xa{B))q By applying Eq. (3) to the composite system A + B and 
then to each system A and B separately and using these definitions, we obtain 
d[{l - q)S*{A)S*{B)] = Y.ayad{X^Yq \ This clearl^leads to quantify the en- 
tropy excess in TT as 5*^^) = (1 - q)S*{A)S*{B) 25]. We see, thus, that TT 
provides also a specific microscopic statistical model for the thermodynamic 
excess entropy and this is the very physical meaning of the nonadditivity prop- 
erty. Each independent system A or i? is analogous to a bulk phase. Despite 
their statistical independency, a thermodynamic coupling exists however be- 
tween both systems. This can be modelled as an interface separating them 
and contributing with an excess entropy 5**-^^ to the properties of the total 
system A+B. K nonextensive system in TT can hence be understood as a two- 
phase-like system at the vicinity of the critical point (since the nonadditivity 
property holds at all scales which means having always significant interfaces 
separating clusters of A from clusters of B). It is interesting that Eq.(9) can 
be recasted as JT" = 5* — 5* — (1 — q)S*S* which is also similar to the Tsallis 
entropy scheme of nonadditivity, see Eq.(2). The SEP J represents the net 



9 



balance between fragmentation or aggregation of the system |l5l.ll6l.ll7l|. This 
SEP is hence to be viewed in TT as a thermodynamic force which can be 
understood formally in terms of an equihbrium of the actual physical nonex- 
tensive system (with physical entropy S*) with its completely fragmentated 
and uncorrelated state in which the system is reduced to a gas of its noninter- 
acting constituents. The gas has entropy —S*. The minus sign comes from the 
fact that this gas represents disintegration, against the cohesion (the positive 
Tsallis physical contribution) which preserves the integrity of the nonexten- 
sive system. This equilibrium is made through an interface contributing with 
— (1 — q)S*S*. Note that Renyi entropy S* is extensive, and hence it describes 
well this completely fragmentated extensive dilute gas (which is an unstable 
and supersaturated one) of the physical system. 

It can be seen also that and S* are related differentially as dj' = — (1 — 
q)S*dS*. This equation is somewhat unusual and has no counterpart in stan- 
dard thermodynamics. If we define / = dJ /dS* we see that f{\A) = Xf{A). 
This means that the chance to aggregate depends linearly on the size of the 
system, which is quite reasonable. 



7 General expression for q and examples 

By using our definition of excess entropy, we obtain 

^ ~ ^ " S*{A)S*{B) " S{A)S{B) 

since S*{A) = S{A) and S*{B) = S{B) through our connection while Eq. 
(6) holds also in NT for the composite system A + B as well as for each 
separate system A and B. In Eq. (10) all the quantities appearing in the 
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r.h.s can be always calculated from NT, which means that there is a way to 
know q from first principles estimations. We now illustrate this assertion by 
calculating q for nonpolar, nonstructured liquids (as argon) at the vicinity 
of the critical point. Liquid clusters develop inside the supersaturated vapor 



phase at all scales. From critical point theory it is well known that the 
surface tension a scales at the vicinity of the critical point (T,. —>■ 0) with the 
reduced temperature = (1 — T/Tc), where Tc is the critical temperature, 
as ~ ctqT/^ (with fi ^ 1.26 in three dimensions, and ctq being the surface 
tension a,nplitud^. The co.elat.on length of phase . (. ^ A «). 5., scales a. 
~ ^oiTj. ^^^^ [2A] where D is the dimension and the amplitude of the 
correlation length. In this limit, the bulk entropies can be written in terms 
of the critical bulk entropic density Sc as Si = ScVi ~ Sc^^, where Vi are the 
volumes occupied by each phase. The excess entropy S^^^ is equal to —ada/dT 



15| where a ~ ^ is the surface area. Replacing all these quantities (taking 
A as the liquid phase) in Eq. (10) and considering D = 3 we obtain 

1-0 = ^ T^'^-i fll) 

where we have introduced the critical wetting parameter Uc = kTd (47r^Q^cro) ~ 



0.78 which is a universal constant for all fluids 



19 



26|. We have shown re- 



cently that uOc is related directly to the dynamics of the critical clusters 



191. 



Concretely, a value ujc = 0.78 means that the critical clusters behave univer- 
sally as a combination of translation and vortex rotational motion jiol. This 
calculation then relates the entropic parameter q to the dynamical quantity 
Uc and the universal exponent fi that can be obtained separately from renor- 



malization group theory |2^. It is to be noted that although renormalization 
group theory is an exact approach given the Hamiltonian, real systems exhibit 
nonideal features due to the complexity of the interactions that can lead to 
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departures from the translationally invariant hamiltonians usually considered 
in these cases. This is particularly dramatic when studying ferromagnetism 

16, 3- 

n our specific example, classical scaling laws are used only here as 
an approximation. It is to be noted that far from the critical point (always 
following the coexistence line) q does not necessarily departs from unity since 
the bulk entropies do not scale with the correlation lengths. Eq. (11) holds 
only approximately as the critical point is approached (as assumed from the 
scaling laws above used). It is to be noted that at the critical point q = 1- This 
is consistent with the fact that the two phases in coexistence merge then into 
a single macroscopic homogeneous one provided that the size of the clusters 
tend to infinity. This means that, in addressing clusters, the thermodynamic 
limit is taken at the critical point (A^ — > oo and V —>■ oo) and hence standard 
thermodynamics should be regained as is the case here. The r.h.s. of this equa- 
tion is vanishingly small for non-structured simple fiuids, as argon, because of 
the large numerical values for Sc- For these fiuids, then, q ^ 1 and this justifies 
the success of the Gibbs method of the dividing surface within standard ther- 
mostatistics in regarding these systems. Departures from this ideal behavior 



2J| in which 



are observed, however, for complex structured liquids as water 
there exists anomalies due to the presence of hydrogen bonding and dipole 
orientation which drammatically complicates the statistical description. The 
presence of structuredness reduce the bulk entropy of the liquid and, further- 
more, interfacial effects strongly increase (as it is known from the value of 
the surface tension of water which surprisingly high when compared to simple 



fluids 



24|). It is then to be expected q to depart signiflcantly from unity in 



liquid water. 



There are many other systems as Ag nanoparticles for which surface func- 



12 



tions can be significantly higher than bulk values 27| which would lead, from 
Eq.(lO), to values of q significantly different from unity. Ferromagnetic ma- 
terials have been also proven to significantly depart from standard thermo- 



statistics 



16l |. and this has been interpreted succesfuUy in terms of a nanoth- 



ermodynamic equilibrium 



16, 3 



Another example of a formal two-phase-like system is a correlated ionic liq- 
uid in the vicinity of a strongly charged macroion 2^. In this situation the 
behavior of the correlated counterions at the vicinity of the macroion is quite 
different to those far from the macroion (constituting a Poisson-Boltzmann 



uncorrelated and disordered liquid) 



28 



29| . For weak coupling, the correlated 



liquid is still a 3D fluid-like one but the rescaled interfacial area available to 
each molecule is affected by correlations and is approximately given by 2S 



29|. It is 



(where S is the coupling parameter entering in the Hamiltonian) 
clear, then, that in the weak coupling regime (0 < S < 1, g < 1) and from Eq. 
(10), 1 — g ~ cS, where c is a small constant (since S''-^-' depends linearly in the 
interfacial area). This is consistent with the previously fitted curve for g(H) 
since in the weak coupling limit we have 1 — g = 1 — 1/(1 + 0. 091S)°'^^ 0.062S 



8 Concluding remarks 

We have established a connection between TT and NT which provides a sound 
and unambiguous physical basis for TT. As a bonus, it has allowed us to jus- 
tify that the y^'s are the correct physical thermodynamic forces (difficulties 
in 3 



pointed out 



are now overcomed) and has led us to two expressions 



that can be used quantitatively to evaluate the degree of nonextensivity in 



13 



Hamiltonian systems. These are Eq. (9), which introduces the connection be- 
tween both formahsms relating thermodynamic smallness to nonextensivity, 
and Eq. (11), which allows for an straightforward evaluation of q from NT. 
A nonextensive system with q ^ 1 can now be understood in the thermody- 
namic limit as one composed of two entangled phases in a critical regime in 
which the (fractal) interface separating them, besides other finite size effects, 
contributes significantly to the total entropy at all scales. 
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